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O ' III these lectures [| we dwell upon the cosmological corner-stones of the Very Early Universe 

^ ' (VEU) theory: Parametric Amplification Effect (PAE) responsible for the generation of Pri- 

'""' . mordial Cosmological Perturbations (PCPs), Chaotic and Stochastic Inflation, Principal Tests 

pq ! of VEU, and others. 
>■ 

a^: 

§ ; 1 Introduction 

11^ . A great success of VEU theory as the theory of the beginning of the Big Bang, is related to its 

Qv^ ! semiclassical nature allowing to operate productively in terms of classical space-time fllled with 

(--| ! quantum physical flelds (including the gravitational perturbations). It (VEU theory) connects 

Q-i| like a bridge the theory of our Universe based on the Friedmann model (FU), with theories 

O ' of Everything (TOE) essentially employing quantized gravity (still very ambiguous). This 

^— > . connection is already realized itself in the important understanding that the quasi-homogeneous 

c^ ! and isotropic state of the Universe on the horizon scale (on one side) and the primordial 

j> I cosmological perturbations which gave a birth to the Universe structure on smaller scales (on 
the other side) are just two features in the low-energy limit of some theory of VEU based on 



Vh ' a model of the inflationary Universe (lU). Up to now we have no alternative to the Inflation 

resulting in the cold remnants which we observe today as the micro and macro worlds, which 
is commonly considered as a necessary element and the probe test of high energy physics and 
any TOE. 

We should begin from the priorities when discussing VEU theory: whether to start from 
cosmological or particle physics standards. The particle physics is not yet flxed well at high 
energies: to follow this direction today means to start from N particle and modifled gravity 
theories where N is a big number, and then to build up N inflationary models based on them. 
For this reason, rather preferable now seems the investigation of cosmologically standard VEU 
theories based of three standard points: 

(i) setting some cosmological postulates within General Relativity (GR); 

(ii) deriving the theory from these postulates, and 

(iii) confronting the theory predictions with observations. 



""^Published in: "Cosmology: The Physics of the Universe", eds B.A.Robson et al., World Scientific, 213 
(1996), and "Cosmology and Gravitation, II", ed. M.Novello, Editions Frontieres, 288 (1996) 
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On this way, which I will follow in these lectures, we are independent of the future particle 
physics and, thus, may try to find some basic properties and principal features of the VEU 
theory (tested by observations!) which stay independent of the particle physics uncertainties as 
well. We know two good examples of such kind theories. The first is the Friedmann-Robertson- 
Wolkel model which has ensured the great success of the observational and theoretical cosmol- 
ogy, just based on the cosmological postulate (the homogeneity and isotropy of the 3-space). 
The second is PAE, the theory of generation of PCPs which gave birth to the Large Scale 
Structure (LSS) of the Universe, just based on the linear perturbation theory in FU. 

Regarding the VEU our goal would be to find the basics for lU in GR which could create the 
Friedmannian region we live in now (which, in its turn, would provide for the initial conditions 
for both our examples, the FU and PAE theories). The answer which we know today is the 
model of Chaotic Infiation based on the assumptions of existence of infiaton (the scalar field 
weakly coupled to all other physical field) and the start-infiation-condition postulating a quasi- 
homogeneous spatial distribution of the initial infiaton in a Compton- wavelength region (i.e., 
in the finite-scale region! A great step in comparison with the FU cosmological postulate !). 

Below, our lectures will be grouped in four parts. The first is devoted to the PAE of scalar 
perturbations in FU (Lukash 1980). Next, we consider the cosmological applications to the 
scattering prohlem. and chaotic infiation (Lukash & Novikov 1992). In the third part we present 
some recent developments of the theory of chaotic and stochastic infiation proposed by Linde 
(1983, 1986). The last Chapter deals with the problem of testing and confronting these theories 
with observations. 

Certainly, we do not pretend to cover all the corresponding references and give a review 
of all recent ideas and speculations in VEU theories. Our main goal here is to present some 
basic properties of VEU which seem today more or less settled and independent of future 
theoretical constructions. We try to consider the simplest mathematic models paying particular 
attention to the physical meaning of the effects considered. Some necessary mathematical 
calculations are given in three Appendices. Hereafter, our units imply c = 8ttG = h = 1 and 
Ho = lOOhkms-'^Mpc-^. 

2 Parametric Amplification Eff^ect 

The formation of the structure of the Universe is one of the fundamental problems of the 
modern Cosmology. The two following properties of the present Universe are very important 
in understanding the physics of its early expansion. 

(i) High degree of homogeneity and isotropy at large scales ( 6p/p^l0^'^ on the cosmological 
horizon) along with a well-developed structure on scales less than 0.01 of the horizon. 

(ii) High specific entropy {N^/Ni, ^ 10^) along with the mean baryon density being less than 
~ 10% of the critic value (l^fe<0.1). 

Property (i) proves that the large scale structure of the Universe (galaxies, clusters and su- 
percluster) stemmed from initially small amplitude perturbations of homogeneous and isotropic 
cosmic medium since it is the small perturbation that may grow up to the order of unity (and 
then form gravitationally bounded object) only when the horizon becomes many times larger 
than its linear scale. We do not know yet whether PCPs formed together with the cosmologi- 
cal Friedmann model at Planckian curvatures of whether they originated in the process of the 



homogeneous and isotropic expansion which is described by the Classical General Relativity 
(CGR). In the first case we have no theory. However, an important point is that the quasi- 
Friedmannian Cauchy-hypersurface is already a classical object after Planckian time. So, if the 
PCPs are made evolutionary and their scales less than the Friedmann-hypersurface scale then 
they are likely to form at a semiclassical stage when the large-scale gravity was governed by the 
CGR equations. We know one example of this kind - it is inflation: galaxy-scale PCPs form at 
the very late stages of the inflationary expansion when the Friedmannian Cauchy-hypersurface 
(which forms the background of our local Universe today) has been already prepared by the 
inflation. Below, we will develop the theory of small potential (scalar) PCPs assuming the 
existence of Friedmann background model. 

Also, it is quite evident that by no process could the Universe be born strictly homogeneous 
and isotropic: there always exist quantum fluctuations of metric and physical fields, the seed 
fluctuations could be of statistical, random character or they might be thermal, etc. An im- 
portant point here is as follows: the inevitable minimum level of seed initial perturbations is 
always maintained by the quantum point-zero fluctuations of a quasi-homogeneous gravitating 
medium which bases the Friedmann spatial slice. 

One of the basic implications of property (ii) is that in the past the Universe was hot and its 
expansion was governed by the gravitational field of intensively interacting relativistic particles. 
At this stage the matter represented to high degree of accuracy a hydrodynamic perfect fluid 
with equation of state p = e/3 (The Hot Friedmann Universe, HFU). 

Small perturbations of such radiationally dominated gravitating fluid are sound waves prop- 
agating through the matter with a constant amplitude (adiabatic decrease of the wave amplitude 
due to the Universe expansion is exactly compensated in this situation by the increase of the 
amplitude gained due to the pressure gradient in the comoving to the wave front reference sys- 
tem). Mathematically, such perturbations in the expanding Universe are governed by the same 
equation as that for the acoustic perturbations of a non-gravitating homogeneous static termal 
bath in the flat Minkowski space-time (conformal invariance). A deep physics is behind it: no 
new phonons, cosmological potential perturbations, can be produced during the expansion of 
the hot Universe. 

For the real generation of PCPs to occur in the early Universe, one has to reject the hot 
equation of state {p = e/3) at some expansion epoch, which can be done in principle only 
in VEU before the primordial heating of cosmic matter to high relativistic temperatures. In 
this case, as we shall see below, the number of phonons is not conserved and new phonons 
(the density perturbations) can be created in the course of the expansion. The mechanism 
is purely classical and works as parametric amplification: energy of large-scale non-stationary 
gravitational potential is pumped to the energy of small-scale perturbations (like new photons 
are created in an electromagnetic resonator when its size changes in non-adiabatic way). This 
effect, which we will generally call "parametric amplification effect", has nothing to do with 
Jeans instability: actually it is the CGR-effect since the typical scales to be amplified are just 
the cosmological horizon (and the effect involves the light velocity and gravitational funda- 
mental constants). Before we present the mathematical formalism the physical meaning of the 
parametric amplification is discussed in the next Section. 

2.1 Physical Meaning of the Parametric AmpHfication 

The theory of small perturbations in the FU was constructed by Lifshitz in 1946. According to 
this theory there are three types of perturbations of the homogeneous isotropic model: density 



(potential) perturbations, vortex perturbations and gravitational waves. We are interested now 
in the first type of perturbations, potential perturbations, because they are related to galaxy 
formation. 

Let us consider the spatially Euclidean background model: 

ds^ = dt^ — a^dx^ = c? [dr^ — dx^\ , (1) 

where a is the scale factor which is a function of time, t and r/ are universal and conformal 
times respectively, 

^ = / - (2) 



a 
It is usual to present the density perturbations bp in a Fourier expansion 

where 5 = 5p/{p + p) (more rigorous definition which makes 5 = 6{t,x) the gauge invariant 
function is given later on). Here, the wavenumber k = \k\ and the physical wavelength A = 
27Ta/k. 

Also we shall use the perturbation scale 

A na .^. 

'' = 2 = T- (*) 

and the horizon scaleQ 

Ih = H-' = ^ (5) 

a 

The latter is a typical scale of the causally connected region during the evolution time scale of 
the a-function. (Dot is the derivative over the universal time (') = ^). 

We consider here only homogeneous and isotropic states of the perturbation fields (ran- 
dom spatial phase fields). Their important characteristic is a power spectrum of the density 
perturbations A^: 

^,{6^,61,) = Al6{k-k'), (6) 

where brackets (...) mean average over the field state, 6{k — k) is the 3-dimensional 5-function 
and (*) is the complex conjugate. This spectrum determines the second correlation of the 
density perturbations 

e(r) = {6{t, mt, X + r)) = JfAf^, (7) 

o 

where r = |r|. All the odd correlations are identically zero. The Gaussian fields which are the 
particular cases of the random fields, can be totally described only by this correlation function 
(all higher order correlations are negligible in linear approximation). Eq. (7) also clarifies the 
physical sense of the power spectrum. The amplitudes A^ are just the corresponding density 
perturbations in the scale interval Ak ~ k around k, they are additive (the total power is a 
sum of the partial ones over all scales) : 



is') = E A 



2 
k 

Akr^k 



^For brevity, we refer to the Hubble scale -ff "'^ as horizon, although it is not technically correct. 



Let us now consider the hot Universe (HFU) when equation of state was p = e/3. The 
solution for A^ is 

Ai = ic,ik)MK)r+ic,ik)MK)r, (8) 



where 



Iff t 1/2 k 



h h v3 

^ , , /sinK cosK — 1\ 
/i K = - cos K + 2 + 

V K K^ ) 

/cosK sinK 
/2(k) = sm k + 2 — 

and Ci^2{k) are the amphtudes of the growing and decaying modes respectively . The sound 
velocity here is /5 = 1/v^ ~ 1, so, k function presents a ratio of the sound pass during the 
cosmological time t to the perturbation scale l^ which is about the horizon to perturbation 
scale ratio. 

The perturbations of the gravitational potential or metric perturbations are as follows (exact 
definitions are given below): 

1-COSK . 

All = Ci tor the growmg mode, 

SlTl /t" 

^2 = C2 — TT- for the decaying mode (9) 

An important feature of these perturbations is the following. Neither growing nor decaying 
modes increase catastrophically in time. Both of them are described by sin and cos functions, 
so, if Ci and C2 are less than unity (and it should be this way, otherwise hi2 would be large at 
smaller time) then both these modes are just sound waves with constant in time amplitudes Ci 
and C2 and with different time phases. 

First of all, this result means that the HFU is absolutely gravitationally stable against small 
perturbations: if initial perturbations are less than unity then they remain small forever till 
equation p = e/3 holds. 

A more elegant proof of this important conclusion may be done with help of the g-scalar 
(Lukash, 1980) which is, generally, gauge invariant combination of matter velocity and grav- 
itational perturbation potentials. Below, we shall see that potential perturbations in FU are 
totally described by this scalar (and back: all matter and metric perturbations can be expressed 
as functions of q) . The physical meaning of the g-scalar easily follows from its definition: for 
large scale, Ik > Ih, Q is mainly the gravitational potential (matter effects are not important), 
while inside horizon, Ik < Ih, gravitational perturbations are negligible and q is just the matter 
(velocity) potential. 

In the HFU the g-field obeys the following equation: 

q + 3Hq-^Aq = 0, (10) 

where H = a/a is the Hubble function and A = d'^/dx'^ is the spatial Laplacian. Transforma- 
tions q = aq, (') = d/dr] = ad/dt reduce eq. (10) to (note, that a ~ 77 for HFU) 

q"-^Aq = 0, (11) 



^In our units the energy density and matter density are the same functions, e — pc^ 



which is just the non-gravitating acoustic wave equation in the flat spacetime (?7,x): 

gfc ~ Ci sin K + C2 cos ft. (12) 

Eqs. (1,10,11,12) indicate the conformal invariance of potential perturbations in HFU and, 
as a result, the conservation of the adiabatic invariant — the total number of phonons, the 
sound wave quanta — which proves the stability of the HFU expansion against small matter 
perturbations Q Note in this connection that gravitational waves have a similar invariance 
property (see Grichshuk 1974) but we do not discuss them here. 

The lumps of the matter in these sound waves start growing only after the equality epoch 
{Ih ~ lO^yrs) when the non-relativistic particles become to dominate in the expansion and 
the pressure falls down in comparison with the total density. This process develops due to the 
Jeans gravitational instability causing the fragmentation of the medium into separate bodies 
at the late stages of the expansion {l^ ~ 10^ — lO^^yrs). 

We shall not discuss here these late processes of galaxy formation. For us the following is 
important: for the formation of large structures (superclusters and clusters), we need a definite 
amplitude of the sound waves ~ 10~^ — 10~^ in the linear scales which encompass the number 
of baryons big enough for these structures formation. So, Ci or/and C2 must be of the order of 
10~^ on these scales. 

It is a very serious demand on the initial perturbations. Indeed, when t is small, k <^ 1, we 
have (see eqs. (9,12)): 

ci<l, 

C2 < K. (13) 

From these expressions we can see that C2 must be extremely small and cannot be of the 
order of 10~^. So, we need in fact the following equations to be met for k, <^ 1: 

(i) Ci > C2, 

(ii) ci ~ 10"^ (14) 

But both of them look very strange. 

Indeed, any general natural initial conditions assume a random time phase state for the 
seed fluctuations 

Ci = C2 < 1. (15) 

E.g., the first eq. in (15) holds for vacuum or thermal fluctuations. More of that, any natural 
fluctuations in hot gravitating medium imply that Ci and C2 are dozers orders of magnitude less 
than 10"^. 

The last point is demonstrated with help of the following example. Let us suppose that the 
origin time of the fluctuations is the Planckian one and let us denote k = 1 for Ipi. Then on the 
galactic scale kgai ~ 10~^^. Now, let us take a thermal fluctuation spectrum at this moment 

^We do not go into further detail about this stabihty effect since we have emphasized it many times in our 
previous lectures. Mention only that the increase in time of the density contrast at k < 1 (see eq. (8)) which 
some people interpret as an instability period, simply corresponds to the period of time of the monotonic change 
of the oscillatory function. (None of the sin-oscillations manage for k < 1). Thus, to speak on the instability 
in this case is as incorrect as to speak on the instability of a mathematical pendulum when it moves, say, out 
of its stable point for the time which is less than the oscillatory period. Returning to our case, note that the 
potential energy of such a pendulum at k < 1 is all in the gravity (see eqs. (9)). The Jeans thinking fails here 
because the k < 1 region is purely relativistic one. 
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with the Planckian temperature and, thus, the maximum at I ^ Ipi. Then the amphtude of the 
perturbations for k < 1 would be proportional to A;^/^ and, on the galactic scale, it would be 
~ 10"'^'^. So, Ci has to be ~ 10~^° and it is 35 orders of magnitude less than we need. 
Our results are the following: 

1) The classical cosmology of the hot VEU has principal difficulties in the explanation of 
origin of the PCPs. Both requirements provided by eqs. (14) for the large scale structure 
formation, cannot be naturally explained within the frameworks of the HFU. 

2) To account for the appearance of the PCPs at the hot Universe expansion period we need, 
as a necessary condition, to reject the p = e/3 equation of state at the VEU stage. The 
modern cosmology provides for a variety of the possibilities of such type: from quantum- 
gravity effects to vacuum phase transitions, cosmic strings, textures, etc. Here we consider 
the most general conditions for the parametric amplification effect appearing in theories 
with one scalar field ip coupled to gravity in the minimal way. 

Parametric amplification means the production of the gravitating potential inhomogeneities 
(PCPs) in a non-stationary gravitational background of the expanding Universe: large scale 
dynamic gravitational field parametrically creates (amplifies) the small scale perturbation fields. 
Mathematically, potential perturbations of FU with a general expansion law are governed by 
the g-scalar which, after the conformal transformations {q -^ q, t -^ rj), meets the following 
equation: 

Opq = Uq, (16) 

where 0/3 = d"^ /drf — /3^A is the light {j3 = 1) or sound type (/3 < 1) d'Alambertian operator in 
the conformal spacetime and U = U{rj) is the effective potential of g-field, which is a function 
of the expansion rate of the FU. Eq. (16) is a type of the parametric equation in mathematical 
analysis capable to amplify the fields with scales k ^ f/^/^ which are usually outside or about 
the horizon size (i.e., in the purely relativistic region). 

Say, for the massless scalar field ip with minimal coupling the effective potential is f/ = a /a, 
thus, the typical frequency is just the horizon one (see eq. (16)): 

"^^^Ml^Hi2-'Jf]'\^H. (17) 

a a \ at J 

For the HFU, {a r^ rj), the effective potential is identically zero U = 0, which reduces eq. 
(16) to eqs. (11,12) considered before. In this case we can define the vacuum state of the g-field 
for all spatial frequencies and introduce, for instance, a standard technics for the scattering 
problem with | in) and | out), vacua, and so on, to see how many phonons are spontaneously 
created during expansion, which are their spectrum, etc. 

Further applications depend on the sign of the second derivative of the initial scale factor. 

The point is that this sign can give us the idea about which scale expands faster: the 
perturbation scale Ik or the horizon Ih (see eqs. (4,5)). Indeed, the first derivative of their ratio 
is just proportional to the second derivative of the scale factor: 

~ a. (18) 

So, if a < at the beginning, then the galactic scales are found initially outside the horizon, 
and the Cauchy initial data should be set up outside the horizon as well. On the contrary, if 
a > 0, then the initial conditions for the scales of interest can be set up inside the horizon. 
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We shall investigate both cases. The qualitative result is as follows: under natural initial 
conditions met by eq. (15) (e.g., | in) vacuum state for the g-field) it is the growing mode 
of perturbations that is finally created in the | out) state due to the parametric effect. So, 
the resulting perturbation field is described by the first line of eq.(14) with the ci{k) spectrum 
depending on the expansion factor behaviour at time period when the parametric amplification 
condition was met {k < U^^"^ ~ aH). 

The Lagrangian theory and the quantization of potential perturbations in FU are considered 
below. The next chapter deals with some cosmological applications. 

2.2 Lagrangian Theory and Perturbations 

Let us consider a scalar field ip = ip{x^) with the Lagrangian density depending on c/? and its 
first derivatives in the following general form: 

L = L{w, if), w"^ = if^iif'' = if^iif^kg'^, (19) 

The action of the gravitating ip field is as follows: 

W\v,g^']= j{L-]^R)^^ d\ (20) 

where Qik and Rik are the metric and Ricci tensors respectively, R = R], g = det{gik)- Variations 
of eq. (20) over (f and g'^^ in extremum give the clasical equations of motions of the ip field 

_y,'M +nz/ = 0, (21) 

and of the gravitational field created by the (/9-field-source 

Ti 

Gik = Tik, Tik = —ipiipk — gikL, (22) 

w 

where n = dL/dw, nv = —dL/dip, Gik = Rik — Rgikf^ and (;) is the covariant derivative in 
metric gik- Note, that eq. (21) can be obtained from the Bianchi identities T^^^ = 0, as well. 

Useful constructions are the comoving (to p-field) energy density and the total pressure of 
the p field 

e = Tikip'^ip'^'/w'^ = nw — L, 

P=lie-T) = L, (23) 

where T = T/. Also, the following equations are valued 

w&^L (9 Inn 



e + p = nw, (3' 



-2 



n dw'^ (9 In w ' 



n w d'^L ^ w d'^L (9 Inn ,„,, 

n op'^ n owoip op 

(Functions m'^ and /3^ can be negative). 

When considering linear perturbation theory p) and g'''' are presented as sums of some known 
functions (the background (°) solution) and small perturbations and hik'- 



Below, we consider the classical backgrounds (eqs. (21,22) are met automatically in (") 
order) and the perturbations can be quantum ones. 

The Lagrangian of perturbation field is got by expanding the integrand of eq. (20) up to 
the second order in and k^^ with the total divergent terms excluded (see Appendix A): 



H/(2)[0,/if]= 1^(2)^/3^^/4^, 



L(2) = L^^\v,ij^) = ^AfW + X'(r' - 1) - 2i;.^^m' + v{y^ - m\ + 2rx)] + 

where v = 0/w, fj = f,j + v{w^i/w), Ui = (p° /w^°\ x = '^iu' — hiku'u^ /2 = 5w/w, ip^ = 
h^ — h6f /2, h = h\ = —ip. (Hereafter all manipulations with indeces are carried out with help 
of the background metric tensors gllJ and g^''^°\ and background index (°) is omitted where 
possible). Obviously, 

^^ X-i^v, -^ = (3-\ + (u + T)v. (27) 



e + p e + p 

The clasical field equations which couple the metric and scalar perturbations, can be ob- 
tained either when the first variations of the action (26) are taken equal to zero or, directly, 
while expanding eqs. (21,22) to the linear order terms. Generally, these equations describe 
three oscillators coupled to each other through the background shear and vorticity (6-order in 
time equation system): one oscillator is the scalar potential perturbations and the other two 
are just two polarizations of the gravitational waves. 0. 

To find the physical degrees of freedom of the perturbation fields and to approach the 
problem of the PGP origin, the following steps have to be developed. 

(i) Gauge invariant functions must be introduced instead of and /ij^ 

The point is that, although the original fields, scalar ip and tensor gik, are genuine by def- 
inition, their decomposition into background and perturbation parts is not unambiguous 
at all. Indeed, if we transform infinitesimally the reference system, 

x' = x' + C, (28) 

where ^* = ^*(x'^) are small arbitrary functions, then the new separation in the coordinates 
£' will take the following form: 

^ = (^(°) (x*) + = (^(°) {x') + 4>, 4> = cp + wi,u\ 

Qikdx'dx^ = Qikdx'dx^, hik = Kk + ^i-k + ^k;u (29) 

where the background metric gi^ has the same functional dependence in the new coordi- 
nates glk (x'') as that in the old ones gl^ {x''). 

To develop the gauge invariant theory one has, first, to expend the perturbation tensor 
hik over the irreducible representations of the background geometry to mark off the scalar 
and gravitational wave polarizations and, second, to find the appropriate gauge invariant 
(i.e., independent of the transformations (29)) linear superpositions of the perturbation 
functions. 



^The vortex perturbations (if any) are standardly found as the first integrals of the 6-order equation system 



(ii) The Lagrangian and Hamiltonian formalisms of the perturbation fields should be devel- 
oped on the basis of the gauge invariant theory. 

An important point here is to obtain the canonical field variables accounting for the 
physical degrees of freedom. 

(iii) Secondary quantization of the perturbations and cosmological applications can be con- 
sidered in connection of the PCP problem. 

Here, we are going to analyze all these points for FU background metrics. In this case 

goi = Ui = 6", Qap = -a^lalS, 

3 a2' 2 a2' 

- + 3H + iy = 0, Vi = V i - (3h{3H + iy)ui, (30) 

n 

e + p 

where •yais = laisix'^) is the metric tensor of the homogeneous isotropic 3-space with 
the spatial curvature K = 0, ±1 (manipulations with the Greek indices are done with 
help of 7q,/3)- All the perturbation types evolve independently of each other in the linear 
approximation since the background shear and vorticity are identically zero. Below, only 
potential perturbations are considered. 

To avoid formal mathematical constructions we try to use here spatially fiat FU (70,^3 = 
Saf^, K = 0, see eq. (1)) and synchronous reference system for its perturbations {hoi = 0) 
if no other cases are pointed out explicitly. The general necessary formulae are given in 
Appendix B. 

2.3 Potential Perturbations in Friedmann Cosmology 

Let us obtain equations of motion of the gauge invariant potential perturbations directly from 
the linearized Einstein eqs. (22), and their Lagrangian from eqs. (26). 

General metric perturbations in the synchronous reference system are presented in terms of 
two gravitational potential, A and B: 

ds^ = dt^ — a^{5ai3 + ha^)dx°'dx^, 

ha/3 = A5a^ + B^a/3- (31) 

Since this metric is governed by the scalar field (19), we have three perturbation potentials 
V, A and B entering the field equations. In fact, only one of them is independent. 

We will consider as the independent one a gauge invariant scalar q which is a linear combi- 
nation of the perturbation potentials. Let us first define the g-scalar and then, using the low 
order equations, relate inversely this scalar to v, A and B. The gauge freedom in the choice of 
these potentials follows from eqs. (29): 



!i,A^A,HF, B^B^FJ% 



V + -, A = A + HF, B = B + fI— + G, (32) 
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where F and G are small arbitrary functions of the space coordinates, and 

r Hi 
2ii = Fui - a\F,k / - + G,k)P^, (33) 

is the most general form of the ^j- vector in the synchronous gauge (Landau & Lifshitz 1967), 
P^ = 6^ — Uiu'^ is the projection tensor. 

Eqs. (32) show that the following function is gauge invariant: 

q = A- 2Hv, (34) 

In fact, this function is independent of any other gauge as well (see Appendix B) which proves 
that g = g(x*) is a 4-scalar in the unperturbed FU. 

To derive the inverse transformations and the equation of motion for g-field we will need 
only the low-order (in time) Einstein equations: 

AA 

6G" = Hh = Se, (35a) 

a^ 

Goa = -A a = (e + p)f ,a, (356) 

SG^a - I^SG:^ = ^{Ci - l^C6^a) = 0, (35c) 

where h = 3A-\-AB and C = A— {Ba^)/a. The first two eqs. (35a,b) are just the conservations 
of energy and momentum, respectively, whereas the last one (35c) states the Pascalian condition 
(the absence of pressure anisotropics for ip-fie\d). In the class of functions under interest eq. 
(35c) have only the trivial solution^: 

C = (36) 

which relates A and B potentials straightforwardly. 

Making use of eq. (35b), we can now express functions v, A and B in terms of the g-scalar: 

^ = ^(Q-|:), A = HQ, B = a-^Q~a-'P, 



Q = 'yqdt, P = a'yqdt, 



= _IL = ^ + P = 3 e + p 
^ H^ 2H^ 2 e ' ^ ' 

Functions Q = Q{x'') and P = P(x*) are the g-integrals over the Friedmannian world line 
dt = Uidx^. They are determined up to the accuracy of some additive functions of the space 
coordinates. This freedom for the Q-function is just the gauge one (see (32)). The P scalar is 
gauge invariant, so its "ambiguity" is physically meaningfull and related to a certain class of 
perturbations of the Bianchi type I model. To prove it we need another relation between q and 
P functions which we are going to get from eq. (35a). 

Let us recover the energy perturbation using eqs. (30, 37): 

^' ^ 3Hv. (38) 



e + p 2Hp^ 



^The formal solution is C = fsiF + gx + e, where /, g and e are functions of time. The last two terms in 
the r.h.s. can be excluded because they do not enter the original hafS functions. The quadratic term can be 
ascribed only to B potential which results in ha/^ ~ fit)5a0, the latter excluded by redefining the scale factor. 
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Now, substituting it into eq.(35a), we have the key equation for g-scalar: 

7a^g = (3^AP (39) 

which is obviously the GR-analog of Poison equation. 

Let us first assume that /9 7^ 0. Comparison of eqs. (37, 39) shows that P(t,x) is specified 
by the g-scalar up to accuracy of additive harmonic function of spatial coordinates P{x): 

AP{x) = 0. (40) 

In the class of the uniformly hmited (in 3-space) functions ha^ ~ P{x)^a/3 the solution is a 
bilinear form with zero trace: 

P{x) = aaf3X'^X^, aaf3 = coust, a^ = 0. (41) 

Thus, potentials v, A and B are reconstructed from the given g-scalar but a partial solution 
that does not vanish under gauge transformations: 

r dt 
t; = A = 0, B = Pix) / — . (42) 

J a'^ 

Appendix C demonstrates that these perturbations (42) are homogeneous and belong to the 
Bianchi type I cosmological model. 

For (3 = 0, function P{x) is arbitrary and eqs. (42) describe the decaying mode of perturba- 
tions. The growing mode is determined by another arbitrary function of the space coordinates, 
q{x). The general solution in this case is 



H ' 

r HAP 

P = q(x) / a-fdt + Pix), 6e = — . (43) 

J a-^ 

So, with all the above said we may conclude that, for (3^0, the g-scalar is totally responsible 
for the evolution of the physical potential perturbations in spatially flat Friedmann model. 
Division over /3^ and differention of eq. (39) (with eq. (37) for the P-scalar taking into account) 
gives the second-order equation for the g-scalar: 

q+(3H + 2^)q-l(Pj Aq = 0, (44) 

where a^ = 7/2/9^. 

Now, we can derive the Lagrangian density for the g-field. Substituting eqs. (37) into eq. 
(26) and leaving out full divergent terms, we have after rather lengthy calculations: 



1^(2) [0,/i^] = W[q] = f L{q)a^dtd^x, 



L{q) = \a' ie - (^) q,.A , (45) 

where L{q) and L*^^) differ each from the other only in the divergent terms. 
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Eqs. (45) evidence that g-scalar is the unique single canonical variable for physical degree of 
freedom of potential perturbations in the FU driven by scalar field of type (19). The Lagrangian 
density depends only on the first derivatives of g-field. However, if a 7^ const then q acquires 
a mass. Endeed, introducing the following transformation 

q = aq, (46) 

we may rewrite the Lagrangian in the form of standard scalar field with square-mass /i^ = —a/a: 

L{q) = l[^- (^) ^^-^'^ - -"'^^'j (47) 

Appendix B confirms eqs. (44,45) for general case. The covariant generalization is as follows: 

(D%);fc = Dik = ^a'{u,Uk + P'P^k), (48) 

L{q) = ^D^'q,q,k, W[q] = j L{q)^/^ d^x, (49) 

where a = (e + p)^/^/2/9if , Pj^ = gik — UiUk, g = det{gik), all Friedmann functions. 

Before going to the next point, we give some relations for density perturbations in different 
systems most frequently used in literature. 

For the synchronous gauge we have from eq. (38): 

^' ' +l{g-HQ). (50) 



e + p 2H(3^ 2 

For the comoving reference system with the synchronized time {vi ~ u. 

5e q 



where 

ds' 



e + p 2Hp^' 


:.|).- 


-a\6ap{l + q)+B^^p)dx''dx^ 


i=t + V, 


x" = x" + y —dt, 


B-- 


q P 



(51) 



(52) 



For the Newtonian gauge (zero shear reference system, B = 0): 



^~' ^ ^3 HP 

+ oW ), (53) 



e + p 2Hf3^ 2 

where 

~2 / HP\ -^ 2 f. HP\ ^ ~a ~(3 

ds = ( 1 ] dt —a [ 1 H I Oapdx dx 

t = t + -a^B, £" + -E'". 
2 ' 2 

For scales in the horizon [kr] ^ 1) all the three expressions for 6e coincide since the leading 

term is the first one {6e ~ g). In the relativistic region [kr] < 1), 6e depends explicitly on 

spatial slice given. 
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2.4 Quantization and Conformal Non- Invar iance 

Taking in mind eqs. (48,49) one can formally consider the g-field as a test scalar field in the 
Friedmann models. It allows for a standard development of the Hamiltonian formalism. 

First, we can construct the Hilbert space of all complex solutions of eq. (48) with the scalar 
product 

(gi,g2) = J J{2dT.i, 

s 

J\^=tD^\qlq2,k-qlu(l2). (54) 

where rfSj is the invariant measure on a Cauchy-hypersurface S. The integral in eq. (54) does 
not depend on S choice because of the 4-fiux conservation law J{2-i = 0. 
Next, the canonically conjugate gauge invariant scalar is introduced: 

a = a{x') = ^IM. = a\ (55) 

oq 

Further steps to the constructing the field Hamiltonian and canonical quantization are as 
simple as that in the case of any other scalar field. We would like to emphasize here two points. 

The quantization is based on the simultaneous commutation relation for the canonically 
operators q and a: 

[q{t, x), a{t, X )] = qa — aq = i\J —g b{x~x\ (56) 

This equation can be compared with the commutator between the velocity potential and 
density perturbation operators of sound waves in the nongravitating static matter, [w, 5e] = 
ib[x — X ) (see Lifshitz & Pitaevski 1978). 

It is worth while rewriting eq. (48) in terms of the conformal coordinates (?], x) for the 
conformal field q: 

Uf^q = Uq, q = aaq, (57) 

where □/3 = ^ — /5^A is the d'Alambertian operator in the Minkowski metric ds'^ = ds^ /a^ = 
drf — dx^. The function U = U{ri) = (aa) /{aa) plays a role of the effective potential of 
scalar perturbations in FU. It is an unambiguous function of the background expansion or, 
more precise, of the scale factor and its time derivatives up to the forth order. Note, that the 
Lagrangian for eq.(57) 

^(g) = ^(g'-/9'g,.r + t/g') (58) 

coincides with L{q) up to full divergent term. 

The total energy of potential perturbations in the Friedmannian space t = const — the field 
Hamiltonian — can be also presented in terms of the conformal field: 

H 



— = a? I Ed^x, 
a J 



E=^,{q' + (3\^r-Uq^), (59) 

where E = E{ri,x) is the local energy density of the g-field. Note, that for the non-gravitating 
matter (or for short wavelengths kr] ^ 1) E is analogous to the sound wave energy density: 

(e + p)^ {f36er 



2{e + py 
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where v = {v'"/a) is the matter velocity. 

Since this formal analogy with sound waves and the fact that PCPs, which are just the 
resulting (after amplification) g-field, are usually found at the beginning of the HFU expansion 
stage, we will call below the g-field quanta as phonons. These cosmological phonons remind 
the standard physics phonons only when phonon wavelength is inside the horizon [kr] ^ 1) [], 
in this case the gravity is negligent and g ~ f (see eq. (34)). For large scales [kr] < 1), matter 
effects are not important and q is mainly the gravitational field potential. 

When scale factor is proportional to the conformal time, g-scalar appears conformally cou- 
pled to FU (see eq. (57)): 

ar^ri/a: U{ri) = 0. (60) 

In all the other cases U ^ and the g-field is conformally non- invariant. It means that 
q interacts with background non-stationary metric, which provides for the spontaneous and 
induced production of phonons in the process of cosmological expansion. 

The secondary quantization of g-scalar results in the following expansion: 



q = J d'kia^q^ + alql), (61) 

where 

a^ and at are the annihilation and creation operators respectively, 

and u = I'kiv) satisfies the following equations 

z/fc + iP'^k'^ -U)Pk = 0, VkT^t - K^'k = i- 

Below, we apply the theory of g-field for VEU. We shall assume | in > vacuum initial 
state for the g-field. To define it explicity, we will consider in the next Chapter two cases for 
1] ^ : a < and a > 0. 

3 Origin of Primordial Cosmological Perturbations 

Here, we consider some cosmological applications of the theory of g-field: the scatterring prob- 
lem and the problrm of the generation of PCPs in chaotic inflation. 

Let us separate explicitely the kinetic and potential terms in the Lagrangian: 

L=p{w)-V{ip). (62) 

From eqs. (23), we have 

Qad = -(e + 3p) = 2V{<^) - {e{w) + 3p{w)) . 



''Rigorously speaking, inside the sound horizon, krj ^ |/3| ^ (see the d'Alambertian in eq. (57)). For 
estimates, we assume in the main text that f3 ^ 1. 

15 



where e{w) = nw —p{w). It is seen, that if V^(v^) > and e{uj) > 0, then the case a < can be 
reahzed only when the potential term is negligible, whereas for the apposite case (a > 0) the 
potential term may play a dominant role. For this reason, we will consider two interesting for 
us asymptotics. 

First, let us suppose that kinetic terms dominates the potential term in general Lagrangian. 
The following theorem can be easily proved in this connection: 

The theory of a real scalar field with Lagrangian depending only on the kinetic term, 

L = p{w), w"^ = (f,iip'\ (63) 

is mathematically equivalent to the theory of potential motions of the ideal fluid with arbitrary 
equation of state 

t \ / \ dpiw) , , ,^,, 

p = p{w), e = e{w) = w — ; p{w). (o4) 

dw 

The 4-velocity of the ideal fluid is a time-like vector {w > 0): 

u' = '±^ = ^^-/w. (65) 

ds 

So, the (f-field acts here as the velocity potential. 

In the other case a > 0, the potential term becomes important and we may decompose 
L-function over small parameter w^: 

2 

L{w, ifi) = -Viip) + —Wi^) + Q{w^). (66a) 

A simple redifinition of (y9-field 

V? ^ / W{ip)dip 



reduces eq. (66a) to the following case with standard kinetic term : 

L = ^^,,^'* - V{^). (666) 

Next Section deals with a general scattering approach for the g-field (the first case can 
be solved only in this approximation). In the last Section of this Chapter we consider the 
Lagrangian (66b) with V{ip) ^ |v^,jV^'*| initial condition. 

3.1 Scattering Problem for g-Field 

First, let us suppose that V{ip) = 0. 

Applying the theorem to the fiat FU, we see that for equation of state p = e/3 the g-field is 
conformally coupled (f/ = 0). Eqs. (59) give the integral of motion 

a r^ rj : H = const, (67) 

which means the conservation of the total number of phonons in the process of the cosmological 
expansion. 

In fact, the phonon numbers conserve at any frequency mode. Let us dwell on it in a bit 
more detail. 
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At relativistic stage (67) phonons are presented by the following choice of functions |: 

1 k 

i/fe = -^== exp{-iujT) , uj = —=, (6^ 

yzu; v3 



where a = const, a/a = t = r] + const, the prime (') is the derivative in conformal time. The 
field Hamiltonian is a sum over all quanta energies: 



Hn^g = J d'kE^N^, (69) 

where Ek = cu/a is the phonon energy and Nj: = ataj: is the operator of the number of phonons 
with physical momentum k/a. The eigenvalues of the mean energy density operator 



are as follows 



^ = 777^ /rf3A:E,(ng+^), (71) 



(27ra)3 J ' y '' ' 2 

where n^ are the occupation numbers of the phonon states. 

Eq.(67) also allows for introduction of the growing and decaying mode operators: 

fuj I a^^ a 
Ci = C,{k) 



^2 ^ C,{k) = \l^ (^^^^] ■ (72) 

In terms of these operators the field expansions have the following form (cf. eq. (12)): 

1 /■ .cr .c^ /^ sintur _ cos cur 



q = ; / d^ke''''^ [ Ci"-^^^^^ + G 



(27r)'^/2 J V UJT LOT 

H = ^ea'Jd'k{\C,\' + \C2\'), (73) 

where \ C \^= CC^ = 0C. 

Now, let us calculate the number of phonons created at a period with some arbitrary ex- 
pansion law governed by the general Lagrangian (19): 

a = a{r]), r]i < r] < r]2. (74) 

We can do it directly taking into account the phonon numbers before [r] < r]i) and after 
[f] > 772) this period and just comparing them (the scattering problem). 

As we have already seen, it is possible to calculate the occupation numbers for any wave- 
length at a linear expansion stage (see eq. (67)). Thus, to solve our problem, we should match 

*We do not go into detail about such standard for any quantized theory things as separation of the Hilbert 
space in the positive and negative frequency subspaces, the Fock space of states, Bogolubov transformations, 
etc. 
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the a-function (and its first derivative a |) by the hnear passes a ~ 77 at the beginning [f] = f]i) 
and at the end (r^ = 772) of the considered period. 

So, in the resulting normahzation (the Planck scale at Planck moment of time is k~^ = 1) 
we have: 

T = 7] = (2t)i/^ V<Vi 

a{r]), rii<V <V2 (75) 

AiT = Ai{r] + r]o), r] > 772, ?7o = const 

where r = T{r]) is the conformal horizon (Hubble) time: 

/,. ^*' 4 /(^ -?)*'■ <-) 

The value of the physical constant 



a 

T= — =7] 

a J Til 



A, = ^(!IZL*) (77) 

a [7] = r]i) 

depends on the average expansion rate at period (74). Locally, Ai-factor can be related to the 
conformal acceleration: 

!L= lim (Ani). (78) 

Let ag and 6^ be the phonon representations (68) diagalizing the Hamiltonian at stages 
rj < rji and r] > r]2 respectively. Then 

h = "feOfc + /^fc'^^p I ttfc 1^ - I A |^= 1, (79) 

where ak and f3k are Bogolubov coefficients. The | in) and | out) vacua are defined accordingly: 

Ofc I in) = 0, bk I out) = 0. (80) 

The Heisenberg state of the g-field is supposed to coincide with the | in) vacuum. It means 
that there are no phonons (i.e., potential perturbations) initially. 

Taking average over the | in) vacuum state, we derive the mean occupation numbers of 
phonons spontaneously created at period rji < r] < ri2, 

{bib^)=nk6{k-k'), nk=\(3k\\ (81) 

and ratio of the energy densities of the perturbations field to the homogeneous cosmological 
field (see eqs. (70,71)): 

The factor Aj"^ takes into account phonon energy cooling during the expansion period (74). 



Calculations of the produced spectrum is also straightforward: 



^We need this condition in order to avoid the unwanted non-physical effects of creation in the matching 
points. For further details about scattering problem, which are standard for any test field theory, see Grib et 
al. (1980), Birrel & Davies (1981), and others. 



qI = (^) (I /3. ? +Re{a,P:e~'^n)- (83) 

The spectrum dependence on the oscillatory exponent means that growing and decaying 
modes are created non-equally (see eq. (73)). Now, we are going to prove that under general 
condition a > {Ai > 1, see eq, (78)) which is most frequently met in the applications, it is 
the growing mode of perturbations that is preferably created by the parametric mechanism. 

Indeed, our initial conditions generally impty {rj < rji, cf. eq.(15)): 

C^ = (|Ci|^) = (|C2r)«l. (84) 

The situation is trivial in case of the acceleration (a > 0) when the perturbation scales 
inflate from inside to outside the horizon. The decaying mode, appearing originally at the 
horizon with the same amplitude as the growing one, decays quickly for larger times while the 
growing mode is frosen. It can be demonstrated with help of general solution of the dynamic 
eq. (46) in large scales (the Laplacian term is negligible): 

q = q-^(x) + q2{x) / — — , (85) 

where ^'1,2(2^) are arbitrary functions of the spatial coordinates. The integral sharply converges 
to a constant in time function 

q = q{x), (86) 

which describes the growing perturbation mode at rj \ Vg |^ q for any expansion law. 

Let us consider in a more detail the case when the initial conditions (84) are set up outside 
the horizon. The general solution (85) helps again. Comparing it with the a-representation 
functions z/^ = i^kiv) which have the form (68) for 1] < rji, and 

z/fc = ^(afce-^-^ + Ae^-0 
for 77 > ?72, we obtain for the Bogolubov coefficients at kri <^ 1 : 

ak = -{A^' + igkA,), pk = -{A^' - iQkAi), (87) 

where g^ = Xi/k — i is the amphfication coefficient, 

— ^ = 1 -2 — / (aa) arj = const. 

\/3 T2^1 J VI 

Obviously, gk = Xi/k for (7^ ^ 1 and gk ~ {uti)~^ for Ai ^ 1. 

Substitution of eqs. (87) into eq. (83) gives the following spectrum for Ai ^ 1: 

k"^ , , sincjT , , 



Qk = 7:— \ gk 



2tt ujt 

The comparison with eqs. (73) reveals easily that only the growing mode is created. 
To clarify the physical meaning of this effect, let us relate directly the Ci,2 operators in a 
and h representations (see eqs. (72,79,87), kri <^ 1): 

cf) = cf) + ^C^\ cf) = Afci^\ (89) 

k 
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So, if one begins with eq. (84) then, in the end, eqs. (89) give for Ai ^ 1 and g^ ':^ 1 

(|CiP)=^^C'2»c2»(|C2p), (90) 

what is just required by the galaxy formation theories (cf. eq.(14)). In fact, the effect is a pure 
game of the mode mixing, it is not the g-field itself that is created but rather the g-momentum 
(the time derivative q). In other words one can say that the parametric amplification effect 
brings about the creation of squeezed state from initially random state. 

So, as we could see, it is not a problem to produce the growing mode of perturbations with 
necessary amplitude. The amplification coefficient is the larger the earlier HFU expansion is 
violated. The typical spectra (88) go like 

Qk = Mk, k < M, (91) 

with the maximum amplitude corresponding to the mass scale M ~ t^^ when the linear 
expansion law was broken for the first time. (Note, that the spectrum (91) decreases to large 
wavelengths in comparison with the Harrison-Zeldovich scale-free spectrum qhz ~ const). For 
k ^ M, the amplification coefficient is exponentially small. 
There are two points concerning eq. (91). 

(i) Initial (vacuum) conditions are set up outside the horizon which requires physical expla- 
nation. 

(ii) To get the expansion factor required for typical scales A; < M to be of the order of the 
galactic scales, it is necessary to ensure the acceleration (infiationary) condition a > at 
period (74). 

In the latter case, the initial vacuum condition must be set up in the adiabatic zone — 
within the horizon — which can be done independently of the expansion law at the beginning. 

3.2 Generation of Perturbations on Inflation 

Let us consider Lagrangian (66) with a potential V = V{ip) G C3 (at least, the first three 
derivatives determined). We shall generally assume it to be a monotonically growing function 
of ip ioT ip > ipo- Without loss of generality we can put (/9o = and V = dV/dip = at v? = 
which, with the symmetric condition V{ip) = V"(— (/?), makes a stable minimum of the potential 
V{ip) at point ip = 0. Under such normalization, our main assumption takes the following form: 

dV , , 

V' > 0, ^ > for (^ > 0. 92 

Also, we define three auxiliary functions of 99 related to the potential derivatives: 

d\a.V c/lnc c/lne 

0? In 93 ' d In y9 ' (i In (^ 

The simplest examples are Vx = m'^ip'^ /2 and V2 = A(/9^/4 where constants m, A are the 
field mass and dimensionless parameter, respectively. Evidently, for the power-law potentials 
Vn ~ v^^", c = 2n = const, and e = / = 0. 
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The background Friedmann quantities and eqs. (30) are 



^' 



2o? 



1, w = n 



3H' = V + 



1 + 



2_r 

02 



-<^, P 



■-|. 



V' 



H = - 



2 ' 



dV 
if + 3HiP + -— = 0. 

dip 



Inflation occurs when 



d = a{H + H^) = -{V - if^) >0. 



(93) 

(94) 
(95) 
(96) 

(97) 



At V > ip"^ the potential energy density contributes doniinantly to the Hubble parameter. 
For our functions (92) it may happen at large (f, and inflationary solution can be got as 
expansion over the inverse powers of | v? |^ 1: 



H 



a = exp 



rH 1 

— (l + -(e-l)cv.-^ + 0(v.-^)), 
(f 3 



Then the inflationary condition can be rewritten in terms of c-function: 
or in terms of potential V: 



dV 



dip 



< Vsv. 



(98) 

(99) 
(100) 
(101) 

(102) 
(103) 



Eqs. (98,99,100,101) are obviously true in a so-called slow-roll approximation which assumes 
friction-dominated equation of motion for y^-fleld (the second- derivative term in eq. (96) is sub- 
dominant). There are two conditions for this approximation following directly from expansions 
(98-101): 



<2, 



e-1 c 



V9- 



< 2. 



(104) 



The flrst condition nearly coincides with eq. (102). From the second condition we have 

d^V 



dip^ 



<2V. 



(105) 
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So, the inflationary solution (98, 99, 100, 101) requires certain analytical properties from 
potential V{ip) (see eqs. (103,105)). We can put it in other terms: eqs. (98, 99, 100, 101) are 
valued for large ip, 

M>^/ = ma^fe /^^]. ((106) 

One can show that, within class of functions (92) satisfying the inequality (106) for 

I 9? |> 9?! = const (107) 

where (fi is a positive root of equation i^i = (/5/(v^i) H ^1^. (98, 99, 100, 101) describe a trap 

(for growing time) separatrix towards which all the other dynamical trajectories |^ with initial 

\ if \> ipi approach rapidly during their dynamical evolution to the stable point ip = ip = 0. 

On the inflationary separatrix (98, 99, 100, 101) ip, H and a vary slowly while the scale factor 

increases nearly exponentially in time, a ~ exp(ift). 

Eqs. (98 - 101) break when the field reaches the point | ^9 |~ yji ~ 1 and the further 

evolution proceeds with damping oscillations around ip = ip = 0. If (PV/dip^ > at (^ = 0, 

then 

2 [2 

H = -r^, ar^t^^'\ (/? = 2W-(mt)-Sin(m(t + to)) < 1, (108) 

o V ^ 

where V ~ vr?'pp' j1 for | </? |< 1, here m and t^ are constants. At this stage the Universe 
expands like a pressureless medium since the average cosmological pressure is exponentially 
small. The medium — coherent oscillations of spatially homogeneous field — is unstable in 
this situation and will decay in particles. The result is reheating and the HFU expansion 
beginning. 

This reheating process although producing some inhomogeneities on the horizon scale ~ fci, 
cannot damage the large scale perturbations created already during the inflationary epoch 
{k< ki). 

To find the postinflationary POP spectrum we must solve Eqs. (61) with parameters (98, 
99, 100, 101) for ip > ipi. 

^'k + (^' - U)iyk = 0, (109) 



aa ' M 2/f2 ' 2a^aH'^ 



2 / 3 



V 



2^1 + ^C3^^2 + 0(^-^)j, (110) 



-'^"For estimates, (pi ~ 1. If there are few roots in eq. ip — ipi{ip) for tp > 1, then the solution (98, 99, 100, 
101) can be broken for some large ip. We do not analyse here such possibility. 

-^^In the phase space {^,<p), initial conditions for the classical trajectories of eqs. (95, 96) are set up on 
the quantum boundary e = 1 which represents a kind of ellipse (circle in case of Vi) around the central point 
p = ip = 0. The radii of this ellipse along ip and (p axes are V~^ [1] and \/2, respectively. 
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where Cijc = c + 6(1 — e), and C3/C = c + 2(e — 1). The conformal time ^ < 0, and initial 
conditions are 

z/fc = — exp(— i/cr]), for k\r]\^ \. (112) 



Eqs. (109,110,111,112) can be solved explicitely by matching two following solutions in the 
overlapping region 

C4V?"^ <k\^\<\ (113) 

where 04/0^ = (c + e — 2)(1 — e) — e/. The first solution assuming the left inequality (113), 
allows for the fZ-potential approximation by f/ = const/r^^ near A; | r] |~ 1 (cf. eq. (110)). Since 
U <^ k^ for A; I ?7 1^ 1, we have for A; | r^ | > C4(/9~^: 

^. = ^y^TM^^^'HA: I ^ I) = ;^(1 - ^) (^"^'' + ^(^"')) ' (114) 

where Hl^\x) is the Hankel function, z/ = (3 + c^Lp~'^)/2. 

Under the right inequality (113), eq.(85) describes the general solution for any Uirf). Since 
the integral in eq. (85) converges sharply in time, we have for /c | 77 |< 1: 

Uk = i7rV2k-'/^aaqk = --^A;"^/^ (^\ g,, (115) 

where constants qk are the PCP spectrum (see eq. (83)). 

Fitting eqs. (114,115) at region (113) gives the following spectrum of the perturbations 
parametrically created outide the horizon: 

where y? = v^^ at the horizon crossing [kr] = — 1) can be expressed directly in terms of the wave 
number 

k = aH = a{^)H{ip). (117) 

The resulting spectrum (116,117) belongs obviously to the growing perturbation mode since 
only for this mode g-field is constant in time outside the horizon. We have already emphasized 
that this property of the growing mode is independent of any expansion law or equation of 
matter state. In particular, spectrum (116,117) holds in large scales for any microphysics 
processes after inflation like phase transitions or reheating. 

Sometimes, people prefer to deal with the power spectrum of density perturbations. Below, 
relation between g-scalar and the coupled density perturbation fleld is obtained in a general 
form. 

Indeed, eq. (44) allows for the growing mode general solution outside horizon: 

q = q{x) + Qit)Aq{x), 

Q{t) = -(^] §(l--/arft). (118) 

So, the comoving density perturbations are (see eq. (52)): 
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Now, it is not a problem to get the desired relation between power spectra: 

We return now to eq. (117) which gives us the connection between k and ip on the inflation- 
ary separatrix (98,99,100,101). This equation can be solved explicitly for a class of so-called 
smooth potentials V{ip). 

Let us introduce smooth potential functions V{(f) for which c{(f) varies even slower than if 

0: 

|e((/?)|<l. (121) 

For such potential eq. (100) yields 

ac:,if-y^exp(-^y ^>^,, (122) 

which shows that ip varies logarithmically in the conformal time (cf. eq. (111)). If Ci = c{ipi) > 
2, then ipi = ci/2 > 1 (see eqs. (106,107)) and 



1/3 



for ip > ifi {\ T] \>\ r]i \) . The substitution to eq. (117) gives for k <ki\ 

Vl = 5 (c, + 41n (^^)) - 2eln(*,/*). (124) 

(the second equality implies k -C ki). 

So, smooth potentials generate the Harrison-Zeldovich types of spectra (see eq. (116)) 
growing only logarithmically to large scales. 

Typical example of smooth potential is a power-law potential 

= 2n " 
which generates the following spectrum: 

g, = (47rn3/2)-ia„y,"+i = 
2 



-(An 



,(n-2)/2^ 



^/2 + ^^'t^^t 



(n+l)/2 



(125) 



where a^ and c = ci = 2n > 2 are constants. 

Important cases are the massive field (n = 1, ai = m): 



...JM4]-,, 



^^Physically, the characteristic scales of smooth potentials are not much shorter than ip. 
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and the A-field (n = 2, 02 



VAJ 


2v^ 


Qk = 


= 




TT 




3/2 

(127) 



For non-smooth V{ip) spectrum q^ can have, in principle, any form depending on given potential 
and the first derivative shapes. Moreover, it is possible to inverse the problem and to find 
potential V{ip) for any given post inflationary PCP spectrum (Hodges & Blumenthal 1989). 
True, some potentials appear to be rather exotic ones, but the result is very important: PCP 
spectra are very sensitive to the potential forms 1^. 

We shall return to the latter problem in the fifth Chapter. But now, let us stress two more 
points in the conclusion. 

Postinflationary perturbations (116,117) are Gaussian with random spatial phases since it 
is the seed point-zero vacuum fluctuations (of the g-field) from which they were parametrically 
created, that are Gaussian by definition. Here, we have no problem with initial conditions for 
the g-scalar because they are determined by microphysics inside the horizon. 

Another interesting point is that most spectra grow with scale growing. It means that there 
exists some critical field (and, thus, the critical scale) for which the corresponding amplitude 
qk ~ 1: 

^ ^ - ^ A;,^ ~ iaH)^r < h- (128) 



n = I 
n = 2 



Say, for potentials (125,126,127) we have 
n>l: ((^i//c)„, ~ 1, fciA,, ~exp(a-2/(-+i)), 

ipcr — 2^77771"^/^, ki/kcr — (m/vr)-^/^ exp(7r/m), 
^cr ^ 3A-l/^ ki/kcr ~ exp(A-^/3)_ 

We shall see in the next Chapter that the Universe on large scales, k < kcr., is globally 
non-linear and it is stochastic (dominated by quantum fluctuations) for ip > ipcr.- 

4 Inflation 

There is no secret that Inflation is a corner stone of the VEU theories. This is not a surprise 
since up to now we have no alternative to the Inflationary Paradigm. 

We are not going to discuss the Paradigm here. There is a lot of reviews and courses devoted 
to the subject (e.g. see Colb & Turner 1989, and references therein). Instead, we would like 
to dwell on the chaotic inflation which, in our view, is the first theory of the kind that can be 
called the cosmologically standard theory. At least, in a sense as this status has the standard 
Friedmann cosmology or the parametric amplification theory. All of them, based on simple 
cosmological postulates which are not directly related to any particular particle physics, can 
explain and predict a lot of obvservational consequences (see the Introduction). 

The goal of inflationary theory is to prepare initial conditions for the standard FU. There 
are the following five items among them. 

(i) Homogeneity and isotropy along with the Euclidean geometry of the spatial slice on scales 
near the contemporary horizon. 



^^Some spectra may violate the slow-roll conditions which makes the inverse problem self-inconsistent. 
"'^'^Physically, non-Harrison-Zeldovich spectra appear from potentials which have characteristic scales less then 
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(ii) The amplitude ~ 10 ^ of PCPs at this shce on galactic to supercluster scales. 

(iii) Reheating sufficient for the primordial entropy production, nucleosyntesis and baryogen- 
esis met in Friedmann cosmology. 

(iv) Small particle numbers (^2 < 1) of the unwanted massive relics created by the Big Bang 
and primordial reheating. 

(v) Small density of the Friedmann vacuum which is the A-term {Q^ < 0.7). 

First inflationary models were rather connected to the specific physical theories and hypothe- 
ses like GUTs, phase transitions, quantum-gravity effects, etc. However, in view of absence of 
the true high energy physics and, which is more important, taking into account a purely cos- 
mological status of the ffist three items above, there was an understanding of the necessity in 
constructing a cosmological standard inflationary theory which could be independent of any 
current speculations about future fundamental physics, on one side, and could account for the 
first three puzzles of FU, on the other side. Certainly, such a theory would not solve the fourth 
an fifth problems which were much more related to the particle physics indeed. 

The ffist theory of such type was proposed by Linde (1983). A basic assumption is that 
potential energy of infiaton ip-Held grows with ip growing (see eq. (92)). The word 'chaotic' 
minds the requirement for a large value of the initial y?- field ((/? ^ 1) which could be realized 
somewhere in spacetime under hypothesis of the chaotic initial conditions. However, we do not 
think that the latter requirement is somehow a problem for the theory at all. 

Below, we dwell on the necessary and sufficient conditions for chaotic inflation and then 
discuss some implications related to the subject. 

4.1 Chaotic Inflation 

Let us dwell on Lagrangian (66) with the potential term of type (92). To start inflation, the 
latter must predominate at y? ^ 1: 

V{ip) >\ ifi.ifi'' I . (129) 

Let us estimate the size of the region where eq. (129) is initially met. 

As we have seen in the previous Chapter the time derivative of initial (/? is not a problem 
regarding the inequality (129) if the spatial homogeneity is postulated, since the inflationary 
solution is a trap separatrix for (/? ^ 1. So, of principal importance is the spatial gradient 
condition following from eq, (129): 

\Vip\<H. (130) 

Eq. (130) can be read as follows: to start inflation one has to prepare a quasi- homogeneous 
distribution of if on scale ~ L = y?/ | V^p \ which is much larger than the horizon scale: 

L>lsi = ifH-^ > Ih. (131) 

The start inflation scale has a physical meaning of the Compton scale of infiaton which 
becomes explicitly clear in case of the massive field {V = Vi): 

Isi = fn~^- 
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Taking it into account for estimates, the start inflation condition (131) is not probably a 
great surprise in a general case as well. 

Eq. (131) deals with the initial distribution of (/9-field. In principle, one can rise a question 
about another (additional to eq. (131)) start inflation condidtion, namely, about initial spatial 
distribution of metric (curvature) on scales less than Isi- However, we will not discuss this 
problem here by the formal reason. From the beginning, we decided to restrict ourselves by 
the case when y^-field is the only source of the metric Qik- The point is that the small-scale 
nonlinear curvature perturbations (if any) assume another source unrelated to the y^-field since 
the latter is homogeneous on scale Isi- So, the curvature born by the i^-field is supposed to be 
homogeneous on scale I si, as well as the y^-field itself. 

Eq. (131) can be interpreted in another way: initial i^-field should be large enough so 
that Ih could be small. However, more stringent constraints for the potential comes from the 
slow-roll conditions (103), (105). Let us explicitly rewrite these conditions in terms of the PCP 
spectrum (106). 

From eq. (103) we have the potential restrictions: 

\/^/2 < lOgfc. (132) 

Eqs (105,116) constrict the spectrum index range: 

dlnqk 



dink 



< 2. (133) 



So, Qk cannot vary from the Harrison- Zeldovich spectrum faster than k^"^ (which is quite 
compatible with the market of galaxy formation theories considered today). 

Eq. (132) puts the direct observational limits on the potential amplitude V{ip) for ip = ip^ 
within the structure scale range k^^ ~ (10 — 10"^)/^,^^ Mpc. Eq. qk ~ 10~^ evidences for the 
weak coupling of (/9-field to the potential V{ip) in this region. 

Endeed, let us demonstrate it when the c-function variation along the scale range can be 
negligent. In this case V = Xip'^, and we obtain a very small value for the coupling parameter 
(in the Planck units): 

A ~ 3c2lO-V~^'^'^ < 10"^ (134) 

Remember that the small coupling parameter is also required for the large size of the FU- 
bubble (to be more than the horizon today, see below). 

Eq. (132) can be used for some other important constraints, e.g., on the reheating temper- 
ature. The radiation energy after reheating cannot actually exceed the inflation energy near 
the end of inflation: 

Prad = l^9*T^H < Vivi) < V{p) < lO'ql (135) 

where v^ = v^fc > v^i {k < ki), g* is the total number of massless degrees of freedom of the 
thermal bath particles. Eq. (135) gives the following upper limit for the reheating temperature: 




(136) 

Making use the microwave quadrupole anisotropy qk ~ AT/T ~ 10~^, we have for the 
standard model {g* ~ 100) : Tfi^H < 10^^ GeV. The latter inequality can be confirmed with 
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help of similar estimate for the gravitational waves produced during inflation (we do not discuss 
this problem here). 

Next important parameter is the Friedmann slice scale Ip, i.e., a typical scale of the part of 
the Universe, created by inflation, which can be approximated by the Friedmann model. Any 
scale in such quasi-homogeneous region is described by eq. (117): 

k = H{ip)a{ip) = H{ip) exp{-N{ip)) [Mpc''^], (137) 

where N{ip) = J Hdt = Nj + Np is the number of e-folds of the Universe expansion from the 
moment when the perturbation was at the inflationary horizon and up to now 

Ni = Ni{ip) = I ^, Np ~ 60. (138) 

After substituting ip = ipcr from eq.(128) we have 

Ip = A;-i = exp ( ^ J ~ exp (x'^) ~ exp(lO^) > 10^^ [cm]. (139) 

The non-linear global Cauchy-Hypersurface which develops in the result of the chaotic 
inflation dynamics, is not built up yet. Nevertheless, we see no principal difficulties to solve 
this problem. The point is that the global spatial Cauchy-Hypersurface cannot exist everywhere 
in spacetime: it breaks in the spacetime regions where the 99-field reaches Planckian densities 
i'-Ppi ~ ^^^^^)j so that the semiclassical approach becomes self-inconsistent. We show in the 
next Section that the latter regions occupy the most part of the physical volume of the Universe 
produced by the chaotic inflation. 

4.2 Stochastic Theory of g-Field 

Let us return again to small scales Ik < Ip where q can be treated as a linear quantum operator 
against the Friedmann background. Equation of motion of g-field is 

g + 3nif g - — Ag = 0. (140) 

where n = 1 + 2a/?>aH, a and H are the classical background functions, see eqs. (93, 94, 95, 
96, 98, 99, 100, 101). As we know from the previous Chapters, the large scale perturbations 
are classical for l^ ^ Ir while the quantum perturbations affect only small scales, l^ "^ Ih- Let 
us separate these two parts of g-field at the inflation period assuming that q is generated by 
quantum perturbation: 

g = $ + F, (141) 

where $ is the classical large scale part of the g-field operator. 

To make this separation explicit let us introduce a notion of the miniuniverse (MU) as a 
part of the actual space-time of the size proportional to the horizon: 

Imu = CH-' > h, (142) 

where ( = const > 1. Evidently, MUs do not expand with the comoving volume. 
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Now, we can define the classical $-field as the mean value of q in MU: 

$ = $(t,f) = f K^{x-x)q{t,x')d^x, (143) 

where K„{r) = (27r)^^/^o"^'^exp(— r^/2o'^) is the Gaussian MU-window, and 

^MU C > 

a = = — - = -riiCv 

a ah 

is the MU-dimension in the comoving x-space, rii = {—ariH)~^. 

Similarly to $ we can define the classical part of the g-field momentum: 

V = V{t,x) = [ K^{x-x)q{t,x)£x. (144) 



The evolution of the coarse grained fields $ and V is governed by the quantum perturbations 
presented by the F-operator: in the inflationary process new and new perturbations created 
inside the horizon inflate, one followed by another, outside the horizon and start contributing 
to the classical fields $ and V when their scales become about (and then larger) than Imu- So, 
F plays a role of the stochastic generator for $, the latter moving like a Brownian particle in 
the gas. More of this, the dynamical equations are similar as well. 

Let us present eqs. (143,144) as the Fourier integrals (see eq. (61)): 



^ = j d?k 0(a,^g. + at^i), 

V= [d'k eia0^, + alql), (145) 



where 

= Q(ka) = e-5'^'"' = f K^irje'^'^d^f. 

Then the original eq. (140) can be rewritten in terms of the classical fields: 

$ - y = /, 



V + 3nHV = ^, (146) 



Hg 

c 

f and g can be called the noise functions (or generators) driven by the quantum fluctuations: 



/ = mHa^ I d'k k^Qia^q^: + atgi), 
g = Ha'J d'k ee{a^,p^ + alp*:), (147) 



where 



Pk = -(Ik + ^r^fc' "^ 



H "'^' aH 

Before we calculate correlators of the noise functions let us introduce the normal modes of 
the classical fields, 

^ = <^ + ^v, v = v, 
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and separate the Eqs. (146): 

$ = /. (148) 

(Eq. (146) for V"-function did not change). Here / = / + vaxgjC^ . The background functions 
n(i) and va^x) are easily derived on the inflationary separatrix ((/? ^ 1): 

n= l + c(l-e)v9^2 + 0(v9-'), 

ni = l-icV-' + 0(^-^), 

m = l-icV' + 0(v^-^), 

mi = l-icV"' + 0(v9~^). (149) 

6 

Eq. (148) coincides with the Langeven equation describing the drift of a Brownian particle 
if the particle coordinate is understood instead of $. 

When comparing windows B and /c^G for the classical and noise functions, we can see 
that in the latter case the main contribution comes from scales / ~ Imv- It is clear: the 
field averaged over the mini-universe, can change its value not before the new perturbation 
reaches the Mf/-scale which happens in a characteristic time (step-time) At ~ Imv- Since the 
perturbations phases appearing on the MtZ-scale are random, the process of the classical field 
change is stochastic. To calculate the characteristic values of this process we must know the 
correlators of the noise functions on the Friedmann hypersurface t = const. 

If the g- field is in the vacuum state then 

(ata^,) = 0, (a^at,) = 6{k - k'), 

and the process is the Gaussian one, so that, the second correlators are quite sufficient to know 
about. In this case the perturbation amplitudes for ip > 1 and ( > 1 are as follows (see eqs. 
(61,115,118)): 

Ik = ^^ (— ) exp{ikx), —q^ = k'^a'^Qk- 

Finally, after the straightforward calculation we have in the main approximation over ip and 

C: 

(//') = (//') = \{{fg) + {gf)) = (gg) = 2D6it -t), (150) 

where D = D{t) = H^{ip/2TxcY is the diffusion coefficient. The (5-function in eq. (150) is used 
instead of each of the following expressions, 

2ij|^ + -J and 12h(^ + ^\ , 

because the halfwidths of the latter bell-functions are about the cosmological horizon (At ~ 
H^^) which is less (by the (^-factor) than the M?7-scale. 

Now, we can introduce the probability distribution P = P{t, $) to find field $ at time t. 
By definition, 

^ Pd^ = 1. 
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Following the standard methods, we can derive the Fokker-Planck equation for this function: 

Obviously, the field dispersion grows in time in this stochastic process. 

Let us take some arbitrary MU at time to with the classical field $o, which we call the 
mother. During infiation the physical spatial volume which belonged to the mother MU at 
t ~ to, expands to larger and larger scales. For t > to, this volume can be covered by other 
MUs (daughters) 1^. The $-field varies from one daughter to another, and the r.m.s. deviation 
(from $o) cr = cr(t) can be calculated making average either by the quantum g-state in one MU 
or over the daughter MUs assembly: 



^2 ^ ^(1, _ 1,^)2^ ^ /■(!, _ |,^)2p^|, 



'i^''--^Jwm^-'- ^'''^ 



It can be confirmed also by the exact solution of eq. (151) 



^"•■^'^Tf^^^pf-^^)- '^^^> 



At the beginning 

a{to) = 0, P{to,^) = 6{^-^o), 

while, during time, the distribution (153) broadened around $o with equal probability for both 
signs of the deviation ($ — $o). The typical one-step-change of the $-field is 




A($) = cr(At) = WV~VC^- (154) 



Before we discuss some implications of this stochastic process, let us consider the necessary 
conditions for the diffusion approach. 

Eqs. (148,150) are generally true if the D-function varies slower than the characteristic 
step-time At ~ CH~^. This requirement is commonly satisfied on the infiationary separatrix. 

The Fokker-Planck approach (151) is less reliable here. Indeed, the notion of the P-function 
assumes that its characteristic change-time should be no larger than the step-time. However, 
regarding eq. (153), it can be marginally so if ( is not too high. Below, we will assume that 

C~i. 

4.3 Non-Linear Inflation 

Mini- universes of size -ff"^(C ~ 1) introduced in Section (4.2), are just the causally connected 
regions of the infiating space. If any two points with constant comoving x-coordinates (i.e., 
expanding with the Universe) belong initially to the same MU then they will manage to ex- 
change the light signals at least once. But if they belong to two different MUs then the light 



-'^^At each step At ~ (^H ^ the mother volume expands by factor N ^ exp(3C), so there are about N daughters 
of the first generation, N'^ of the second, and so on so forth, inside the volume. 
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signal sent from one point will never get the other one. An important consequence is that each 
MU expands in time independently of any other. 

It means that any MU can be chosen as mother regarding the next generations of the daugh- 
ter MUs it produces. In its turn, any daughter taken at a moment t ^ to, although created by 
some mother at t < to, is the mother itself for t > to- This picture has neither beginning no 
end. Actually, this boiler of MUs is eternally self- reproducing inflationary Universe. 

Let us consider some physical volume expanding with the comoving space. New and new 
MUs are created inside the volume during the evolution. We can connect by the time-like 
tracks causally related MUs (mother-daughter, mother-daughter, etc.). We saw in the previous 
Section that there is an equal probability for both signs of the field deviation A<l> to be found 
along any track from the past to future. Since any MU develops independently of the previous 
history and its neighbours, we can forget about the seed mother field ip^°'' and try to find the 
current classical quasi- Friedmannian field <fMU driven the given MU on the track and applicable 
only to this MU. Certainly, ipMU is the local (y9-field renormalized each time by the classical 
part of g-field (see eq. (141)). 

Technically, we can use Newtonian gauge to find ipMU, since this particular frame most 
closely imitates the local Friedmannian expansion (see eqs. (53)): 

dslju = (l - —^^u) ^^' - «' (l + —Pmu) dx\ (155) 



where 



The local time is 

H 



Pmu = / a'j^dt. 

/TT 
—PMudt. 



and the local expansion factor is 



TT 

a-MU = a ( 1 + —Pmu 



In both equations we disregarded the dependence of the Pmu function on the ^-coordinates 
within MU. In doing so, we can easily recover the (^Mc/-field from background eqs. (95). In 
the main approximation over y? ^ 1, we have: 

^"^ 1 + -(i<I>--Pmu). (156) 



(^(°) c V a 

Below, we give only qualitative ideas about some results of this investigation. 

The most interesting question arisen is as to which densities of the y^-field most tracks lead 
during stochastic inflation? The answer is a production of two factors: the probability to find 
a certain y9-field on one track, and the number of tracks earring given (/9-field. As far as the 
first factor is concerned, the classical monotonic decrease of (f when sliding down the potential 
V{ip), is the smaller the larger the field is: Aip^i ~ ~c/'^ for one step-time At ~ H^^. On 
the other hand, the field stochastic change due to the quantum perturbations proceeds in both 
directions of (/? with amplitude | Acpst |~ H growing to higher field values. So, for large enough 
(p{(p > (c^/A)^/'-^"'"'^^) the quantum stochastic process predominates, so that, at each step half of 
the created daughters have higher (f than their mother. But the total number of tracks created 
per unit time grows to higher fields as well, dN/dt ~ 20H. So, the majority of tracks leads 
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to high densities, thus, the largest part of global physical volume is occupied by the Planckian 
field density, i.e., by the space-time foam. 

We can only guess what is happenning there, in these most typical high density states of 
the global Universe, — the notions of the space-time and infiaton break, mutable transitions to 
different physics, signatures, dimensions and other conceivable and inconceivable worlds may 
occur. In fact, we can only say that the inflation states with densities less than the Planckian 
one, are non-typical and very unprobable ones in this really chaotic Universe dominated by 
the sea of quantum fluctuations. The stochastic regime considered above is just a part, the 
semiclassical part of this sea, where the space-time is already classical while the inflaton is still 
dominated by quantum fluctuations. Evidently, such regions decouple occasionally from the 
space-time foam and exist independently during some period of the classical time. 

A very important conclusion is that inside these semi-classical regions there exist some very 
few tracks which lead occasionally (through the random stochastic process) to lower and lower 
densities of 99-fleld. When the latter becomes below the critical quantity (y? < (c^ /XY^^'^^'^^), 
quantum fluctuations are not able any more to increase (f in the created daughter- MUs 
and the successive inflation continues with a monotonic decrease of the inflaton. Now, the 
quantum fluctuations are responsible only for small density perturbations varing slightly from 
one daughter to other. So, the result will be the adiabatic Gaussian perturbations (with the 
amplitude decreasing to smaller scales) against the Friedmannian background patch surrounded 
by the non-linear chaotic Universe. 

Let us emphasize two other points in the conclusion. 

The global Cauchy-Hypersurface does not rigorously exist. It can be constructed near the 
Friedmannian patch in a space-time region restricted by the Planckian densities. In fact, such 
a non-linear solution describes just a temporal island formed with a very small probability in 
the chaotic space-time foam of the Universe and suitable for life. In this connection, we can 
mention that we do not think it is worth while putting seriously the question as to how general 
are the start-inflation conditions for the chaotic or other inflation theories from the point of 
view of general solutions of the CGR equations? In our opinion, it is quite enough that the 
probability for creating a low-density world where life can appear, is non-zero. We have no 
time to discuss this subject in more detail here. 

The next important point is as follows. There are the inflaton quantum fluctuations in 
the chaotic inflation theory that are the reason for both, the non-linear global structure of 
the Universe as a whole and the adiabatic density perturbations responsible for the large scale 
structure formation in the Friedmannian patch of such a Universe. Therefore, we can test the 
inflationary theories just investigating the spectrum of PCPs and then reconstructing the global 
structure within the theory frameworks. This is, probably, the only informative cosmological 
channel to learn anything about the features beyond optical horizon, as well as the fundamental 
physics parameters beyond direct experiment. We can see that the problem of testing inflation 
becomes part and parcel of the VEU theories. We are going to discuss it briefly in the next 
Chapter. 

5 Testing VEU 

There are few important cosmological predictions coming from the very early inflationary epoch. 

Among them are the total energy density in the Universe and the spectrum of adiabatic PCPs. 

The former quantity is to be equal to critical density up to accuracy of the PCP amplitude 



33 



on the contemporary horizon: 

f]toi = l±5H, \K\<6h. (157) 

This amphtude Sh = Sk{k = Ho) can be easily estimated by the quadrupole anisotropy of 
the microwave background radiation: 

6h ~ 10~^ 

The real dynamical density is close to the critical one with the accuracy ~ 30%. However, 
even if eq. (157) were confirmed with a much higher accurary by future observations, it could 
not, unfortunately, tell us much information about the inflation principal parameters. In this 
respect, more informative is the PCP spectrum generated by last stages of the inflationary 
epoch before the beginning of FU. 

As far as the chaotic inflation is concerned, the postinflationary PCP spectrum is very 
sensitive to the potential form (see eqs. (116, 117). On the other hand, the shape of the 
potential energy of the inflaton found in a given Particle Physics, must be unambiguously 
fixeds by the fundamental Lagrangian regarding all the particle fields and interactions. 

Thus, a principle test of inflation is the large scale structure of the Universe. Its analysis 
allows for restoring of the postrecombination PCP spectrum on galactic to horizon scales. This 
part of the spectrum is obviously related to the postinflationary PCPs which, in their turn, 
depend directly on the inflaton potential within the field interval responsible for the scales 
mentioned above (see eq. (117)). 

Two following questions arise from this consideration: 

* How to relate the postinflationary and postrecombination spectra? 

** Which are the basic cosmological observations now pouring the light on the PCP spec- 
trum? 

Below, both topics are discussed very briefly. 

5.1 Transfer Functions 

The point is that any real confrontation of theory with observations can be done only within 
the framework of some cosmological model allowing to transfer PCPs from times when they 
were produced and up to the moment when they entered the non-linear evolution to form the 
hierarchy of the objects observed. The principal parameters of the model are those of the dark 
matter components running out the gravitational evolution of PCPs. 

Regarding the gravitational impact to PCPs there are two components of the Dark Matter 
Cold (CDM) and Hot (HDM) that are important. The cold particles may be heavy relics 
{mc ^ lOeV^) or coherent axions which behave like a non-relativistic medium. The hot particles 
are those like massive neutrino with the equilibrium particle density and according restmass 
{niH < 10 eV). These two components evolving very differently in the past are both non- 
relativistic and maintain the critical density now: 

nc + nH + nA = ntot = i- (i58) 

Baryons can be added in Qcj the third term which is the energy of physical vacuum cannot be 
totally excluded today < Q\ < 0.7. 
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The rest dark matter, which we call as the z/-component, does not contribute crucially to 
eq.(158) and consists of relativistic and semirelativistic weakly interacting particles {rrii, <^ 
leV). We can characterize //-particles by there total number N,^ with respect to the relic 
photons: 

For the standard CDM and HDM models u = 0.4 and u = 0.3 respectively, counting three 
or two sorts of the massless neutrino. Generally, //-particles include gravitions, light SUSY and 
other hypothetical ...inos probably existing in the Universe. 

So, in the simplest matter dominated case {fltot = 1, A = 0, stable particles) we have two 
free parameters, Qh and u, both ranging from zero to one, which determine the past history of 
PCPs beginning from Inflation. The goal is to find the ratio of the final to initial PGP spectra 
as a function of these two (or more in a general case) parameters. This ratio is called the 

transfer function, 

„{/) 
Tik) = \, (160) 

where postinflationary spectrum q^ coincides with the function q^ from eqs. (116, 117), g^ is 
the postrecombination PGP spectrum responsible for LSS formation in the Universe. Evidently, 
T{k) does not depend on the inflationary period, it is a functional of the Friedmann model 
dynamics from the beginning to our days. 

T{k) is equal to unity for very large scales (T(0) = 1) and then decreases monotonically with 
k growing. It still remains to be about unity up to some characteristic scale leg which coincides 
with the horizon at period at equality of all relativistic and all non-relativistic component 
densities. The further T(A;)-fall-shape to shorter wavelengths is an intrinsic property of the 
model (we refer this subject to the special courses and lectures). 

For the standard models leq ~ 30h^^ Mpc but, for arbitrary u, l^q grows with v growing like 

As we see, the resulting spectrum qi is a sensitive function of both fundamentals of the 
early Universe, the postinflation PGPs and dark matter composition. So, the investigation of 
its direct creature — the large scale structure of the Universe — cannot be overestimated today. 
Another principle test — a laboratory detection of the dark matter particles — is not discussed 
here. 

5.2 Observations 

Since there are few special co ureses devoted to this subject we only briefly outline the hot spots 
of the confrontation between observations and theory important for us. We see tremendous 
importance for the modern cosmology of two groups of experiment nowdays: AT/T observations 
both on large and small scales (0^1°) and direct observations of the distribution and evolution 
of the hierarchies of LSS. 

The point is that both experiments confront and complement each other. 

If AT/T upper limits and detections which become available now, make us to lower down 
the primordial perturbation amplitudes on scales larger than / > 10h~^ Mpc, then LSS needs 
for its existence high enough cosmological perturbation amplitudes on scales / ~ 10 — 100h~^ 
Mpc. For the most theories of galaxy formation the gap between these two requirements is 
quite negligible. Say, within the Gaussian perturbation theories any reasonable asumption for 
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large superclusters and voids to be more or less standard phenomenon in the visible Universe, 
leads inevitably to the AT/T prediction levels on degrees of arc capable for current detection. 
It brings a very great optimism to obtain large scale primordial perturbation spectrum directly 
from the observations with a high degree of accuracy. 

The current situation with AT/T is well known. For our case of Gaussian PCPs we may 
directly relate our scalar q on the last scattering surface at recombination with the map of the 
temperature anisotropies on the selestial sphere n = n{9,ip). Endeed, let us decompose the 

latter in spherical functions: 

AT 

-7fr{n) = Y^ aimYi„,{9, v?). (161) 

Then, after simple calculations, the temperature correlation function takes the following 

form: 

AT AT 
C{a) ^ { — (n,) — in,)) = J^QPiicosa), (162) 

where nin2 = cos a, Pi{cosa) are the Legandre polinomials, (...) is the average over the field 
state, and 

Generally, there are three main sources of the primordial temperature anisotropies: due to 
fluctuations of the gravitational potential, matter density and velocity perturbations. For the 
vanishing pressure they may be reduced to the following expression taken at recombination: 

AT 1 1 

-jrin) = -q{xr) + -V - nvr. (164) 

The first term (Sachs- Wolf effect) dominates for a > 1°(/ > 100 h^^ Mpc) which yields the 
direct connection between C{a) and the PGP spectrum g^: 



1 1 fllJ^k), C{a) = ^l ^qlUm, (165) 



where J; are the Bessel functions, /? = 2 sin |. 

The most important is the GOBE detection for a ~ 10° evidencing the consistency with 
the HZ-spectrum on very large scales, / ~ 1000 h^^ Mpc: 

Alr^k^qlr^k"", a = 4±0.2, (166) 

where (5^ = (1 + z)"^ / ^A| is the mean square perturbation of density. 

As for the LSS data, we have many independent indications today for existence of the large 
scale structures up to a typical scale Ils — 100 h~^Mpc. The most important data come from 
the Great Attractor {z < 0.03), pencil beam galactic surveys {z < 0.3), existence of large groups 
of quasars {z G (0.5,2)), and spatial distributions of clusters of galaxies {z < 0.1). These LSS 
experiments indicate the following estimate for the spectrum of Gaussian density perturbations 
within the scales / G (10, 100) h~^ Mpc: 

A^~F, 7 = 0.9 ±0.2. (167) 

The consistency with eq.(166) gives us the obvious evidence for the presence of a real feature in 
the power spectrum at the supercluster scale ~ 100 h~^ Mpc; it may a change in the spectrum 

36 



slope from HZ (at / > 100 h~^ Mpc) to the flat one (/ > 100 hMpc). This "signature of the 
God" requires its explanation in physics of the very early Universe. 

We do not discuss here modern status of the cosmological model and observational tests 
refering to special reviews (e.g. |astro-ph/980321^ ). 



6 Conclusion 

As we have seen, still in the absence of high energy physics, we may successfully develop the 
theory of VEU on purely cosmological grounds and come to important conclusions about the 
spectrum of PCPs capable of current testing by observations. On the other hand, we are 
very close to recover the postrecombination PGP spectrum directly from the observations, 
both AT/T and LSS (especially its evolutionary aspects), and thus to reconstruct the true 
cosmological model and make the exciting link to VEU physics. Both confronting branches — 
theory and observations - develop fruitfully and make us hope to find the principal answers on 
the evolutionary model of the Universe during our days. 
This work was partly supported by INTAS 97-1192. 

Appendix A 

Here we obtain the Lagrangian L^"^^ (see eqs. (26)). 

Asuming that eqs. (25) are exact ones, all the auxiliary quantities are expanded to the 
second order in and /i*'^: 

6L = n'-^'Udw - z/(°)0 + r— + —{{—)^ - m^^) = 

w 2w (3 

= n^°'^w^''\x - ^^°^v + ^Xv + \{I3'^ - 1) + \viv' - UiVkh'^ - ^-m^v") 



^\kul 



Sgik = hik + hiih^, \n\-j-^\ = h + -hi hj^, 

1,.., 1 



.9 



SV^=^i-9^"^ {h + -hthi + -h'), 

where Sf = f — f^°\ v = (f)/w^°\ vi = (l)^i/w^°\ P'^ are the Ghristoffel symbols. 
The substitution to eq. (20) yields: 



(^ - ^^) Ji) = ^'"' - W"' + ^(1 + \f^^^- 



gio) 
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L(2) = -nw{viv' + x^(/5"^ - 1) - 2t;i^^M^ + uv^Ij- 



4^ 2 



-m^t;^ + 2rt;x) + -(L° - -i?(''))(^,,^^'= - -^2^ + 



1 

2* 



If background metric satisfies the classical equations, 

then the substitution to L*^^) gives eq. (26). Note, that the linear terms in and h^^ (see eq. 
(20)) prove to be zero since the background equations are met. 

Appendix B 

Here we derive general relations between g-scalar and the perturbations in arbitrary reference 
frame. 

Let us decompose perturbations (25) over the irreducible representations of potential type 
in general Friedmann model with the metric gik and 4-velocity -u*: 

T; = X + i(C + DX), 

hik = Ycik + Zgik + (CM(i);fc) + D-ik, (169) 

where 4-tensor Cik = 2uiUk — gik has Euclidean signature and subbrackets mean the sym- 
metrization. Functions X, Y, Z, C, D are coefficients of the linear decomposition. X, Y and Z 
are gauge invariant 4-scalars. C and D are arbitrary functions specifying the gauge feedom of 
perturbations in eqs. (25) (see eqs. (28, 29)): 



The Einstein linear equations can be decomposed as a 4-tensor over the irreducible representa- 
tions as well: 

h^L,i - 2/^(i;fc)/ + h,ik + {p- e)hik+ 

+ (e + p){x{P~^ - l)eifc + 4:V(iUk) - 2uvgik) = 
= Etik + Fgik + (/M(i);fe + J-ik = 0, 



where 



E = Y^ - 4HYiu^ + {p-e + 4{H + H^))Y+ 

+ (e + p)(x(r'-l) + 2(4iJ + z/)X), 

F = Z:l + {p- e)Z - AH{Y^iu' + HY - {e + p)X), 
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Ui-k = H Pik = -H {gik - Cik) , H = -u\, 



I = -2(Y^iu' + HY - {e + p)X), 
J = 2Z, 
and X = w^^{wX)^iv} — (F + Z)/2. Here, the auxiliary relations 

-H{gik-e^k), ^=3' 

and some other background formulae were used. 
From J = we have 0: 

Z = (170) 

From I = F = the following relation between X and Y scalars is obtained: 

Y + HY = {e + p)X. (171) 

For the spatially flat model {K = 0) g-scalar is given by the linear superpositions of the 
gauge invariant functions: 

q = -{Y + 2HX) (172) 

Obviously, it is the real 4-scalar (as well as X and Y) independent of any reference frame. 
The inverse transformations follow from eqs. (171, 172): 

X^i(f-I). Y^-Hp, (173) 

2 \a H / a 

where P = J a'-yqdt, 7 = —H/H"^. 

From E' = we have the key equation |^ 

a37/5-2g = AP. (174) 

The rest is to prescribe the potentials C and D for different gauges. 

Projecting eqs. (169) on the Friedmann reference system we have for an arbitray gauge: 

v = X + -F, hoo = Y + F, 

Ka = \^,a. hi = Abi + E'f , (175) 

where F = C + I), ^\) = a^ipja^) + F, A = HF - Y, B = ~D/a^. So, the most general 
definition of g-scalar in terms of 3-potentials is the following (cf. eq. (34)): 

q = A- 2Hv. (176) 

For the orthogonal gauge {ip = 0): 

C = 2a{aB), D = -a^B, F = a^B. (177) 

The next examples specify the function B in eqs. (177). 
For the synchronous gauge {hoo = ip = 0): 

B = a-^Q-a-^P, Q = f -fqdt. (178a) 

For the comoving gauge {v = ip = 0): 



1 
Jc 

For the Newtonian gauge B = ijj = (cf. eqs. (51, 52, 53)) 



B = -^- a-^P (1786) 



-'^^The scalar J is generally connected to the anisotropic pressure which is zero for (/3-field 
^'^Functions (e + p)x and (e + p)X should be excluded from E with help of eq. (171). 
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Appendix C 

One can easily verify that eqs. (42) are just the hnear expansion terms of the exact solution 

dt 



ds' 



df - a^ exp{2Taap)dx''dx^, t = f —, (179) 



where the function a = a{t) can be found as follows: 



1 1 9 



For t ^ 0, we have the Kazner asymptotic: 



a' 



3 



At, (7„^ ~ diag(t2p<^), (180) 



where Kazner exponents (pi + P2 + P3 = Pi + P2 + Ps = ^) ^^^ obviously related to the eigen 
values of matrix a^^: 

det{aa/3 - XSafs) =0, Xa = Aipa- -j ■ 

Eqs. (170) describe solution for the Bianchi type I model with comoving space. The infinite 
scale vortex and gravitational- wave perturbations lead also to eqs. (179). Note, that the 
infinite scale perturbations although causing the expansion anisotropy (shear), do not perturb 
the spatial curvature and density perturbations {6e = Ua = 0) ■ It happens only in spatially flat 
Friedmann models. 
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